A model for the formation energies of alanates and boranates 
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We develop a simple model for the formation energies (FEs) of alkali and alkaline earth alanates 
and boranates, based upon ionic bonding between metal cations and AlHJ or BHJ anions. The 
FEs agree well with values obtained from first principles calculations and with experimental FEs. 
The model shows that details of the crystal structure are relatively unimportant. The small size 
of the BH4 anion causes a strong bonding in the crystal, which makes boranates more stable than 
alanates. Smaller alkali or alkaline earth cations do not give an increased FE. They involve a larger 
ionization potential that compensates for the increased crystal bonding. 



The large scale utilization of hydrogen as a fuel cru- 
cially depends on the development of compact storage 
materials with a high mass content of hydrogeufi Over 
the last decade alanates and boranates have been stud- 
ied extensively because of their potential use as hydrogen 
storage materialsfi^ These materials consist of a lattice 
of metal cations and AIH4 or BH4 anions, respectively. 
The ideal hydrogen storage material should have a high 
gravimetric hydrogen density, which requires the use of 
light metals. Moreover, the formation energy (FE) of 
such a material has to be such that it is stable at room 
temperature, yet it has to decompose at low tempera- 
ture to release its hydrogen. In principle a large variety 
of alanates and boranates can be synthesized by chang- 
ing the metal cations, which can be used to tune the 
formation energy.'* 

Since synthesis is a very time consuming effort, there 
is a need for a materials specific theory with a pre- 
dictive power for the FE. At present the state of the 
art is formed by first principles calculations based upon 
density functional theory (DFT). Several papers have 
been dedicated to trends in the DFT FEs of alanates 
and boranatesi^i^i^ii There exists a surprising variety 
of crystal structures among these compounds. In DFT 
calculations the crystal structure with the lowest energy 
has to be searched for each compound, and the cell pa- 
rameters and the atomic positions have to be optimized. 
This procedure also makes DFT calculations a very time 
consuming effort. A simple theory would help to under- 
stand the trend in the FEs of alanates and boranates. 

Our aim is to construct a simple model for the FEs at 
K of alkali alanates and boranates (MAH4, M = Li, Na, 
K; A = Al, B) and of their alkaline earth counterparts 
(M'(AH4)2, M = Mg, Ca), avoiding the use of the ac- 
tual crystal structure. We assume that these compounds 
can be described by ionic bonding between M+ or M'^+ 
cations and AII4 anions. Our model for the FE, Ai?/, is 
based upon a Born-Haber cycle,— 



Starting from bulk elemental solids and H2 molecules, 
-E'eicm is the energy required to atomize the solids and 
the molecules. The £^eiem of the bulk solids are listed in 
Table HI We use a value of 4.48 eV for the dissociation 
energy of H2 ^ 

The second step is to create M+, M'^+ and AH^ ions 
from the atoms, represented by the energy -Eions- The 
contribution to E'ions from the M+ ions is simply the first 
ionization potential (IP) and from the M'^+ ions it is 
the sum of the first and second IPs. The numbers S/P 
are given in Table |T1 We calculate the contribution to 
-Eions from the AH^ anions as follows. First an electron 
is added to an Al or B atom, which lowers the energy 
by the atomic electron affinity (EA). These atoms then 
have four valence electrons that are used to form covalent 
bonds with four hydrogen atoms. Using EAs of 0.44 eV 
and 0.28 eV for Al and B, and 2.91 eV and 3.45 eV for 
the Al-H and B-H bond strengths,^ we calculate FEs of 
-12.08 eV and -14.08 eV for the AIH4 and BH4: anions. 

The final step consists of constructing the crystal from 
the M+ (or M'^+) and AIl4^ ions, which is represented 
by the energy E'crys- We use a simple Born model for 
the potential between cations and anions. It consists of 
an attractive Coulomb potential between point charges 
at the centers of the ions plus a repulsive short-range 
potential oc r~", where n is the average Born exponent, 
i^crys is then given byi^ 
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where Za = — 1 and Zc = +1, +2 are the valencies of the 
anions and the cations, respectively, and tq is the short- 
est cation-anion distance in the lattice. Mc represents 
the Madelung constant, which depends upon the type of 
latticeJ^ 

Note that in all compounds considered here the AIII4 
and BII4 anions have a tetrahedral geometry. However, 
in the Born model of Eq. ^ , we have approximated these 
tetrahedra by spheres. Our motivation for this is that we 
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TABLE L Ionic radii, Rion (A), summed ionization poten- 
tials, E/P (eV), dissociation energies for the elemental bulks, 
EdiB (eV), and the Born exponents, n. Values are taken from 
Refs.iand [U 





Li 


Na 


K 


Mg 


Ca 


B 


Al 




0.90 


1.16 


1.52 


0.86 


1.14 






EIP 


5.39 


5.14 


4.34 


22.67 


17.98 






-Edia 


1.64 


1.08 


0.93 


1.48 


1.81 


5.81 


3.38 


n 


5 


7 


9 


7 


9 


7 


9 



are interested in a simple model of -Ecrys without having 
to take into account the full details of the crystal struc- 
ture. The cation-anion distance then is the sum of the 
ionic radii of the cation and the anion, ro = rc -\- rA- 
Since the cations we consider are mostly octahedrally co- 
ordinated, we use standard ionic radii rc of 6-fold co- 
ordinated alkaline and alkaline earth ions, see Table [Bii 
As the radius of the anions ta we use the Al-H and B- 
H bond lengths, which are 1.62 A and 1.20 A, respec- 
tively. rA then roughly corresponds to the average of the 
maximum radius of an AH^ ion (A — Al, B) and the 
minimum radius, which is the radius of the central atom 
A. The values obtained for the cation-anion distance tq 
then correspond to the average of the cation- Al/B and 
cation-H distances in the crystal. 

Avoiding the full details of the crystal structure also 
leads to using average values for the Madelung constants 
Mc in Eq. The alkali alanates and boranates have 
an AB type lattice, where A is the alkali cation, and 
B is the boranate or alanate anion. The variation of the 
Madelung constant over different AB lattices is relatively 
small, so we use an average value Mc — 1.76. The root 
mean square deviation (rms) averaged over all AB lattice 
types is 4%. A similar reasoning holds for the alkaline 
earth alanates and boranates. They have an AB2 lattice, 
whose average Madelung constant is Mc = 2.40 with a 
rms deviation of 4%. 

Fig. [T] shows the most important results, i.e. the FE 
calculated with the model represented by Eqs. H]) and 
compared to experimental valueSfi^ For some of the 
materials the experimental FE is not known. There- 
fore we have also performed first principles DFT calcu- 
lations. We use the projector augmented wave (PAW) 
methodfi^ii^ and the PW91 generalized gradient approx- 
imation (GGA)fi^ as implemented in the Vienna ab ini- 
tio simulation package (VASP) ]^^i^^'^^i^^ To integrate the 
Brillouin zone we apply the tetrahedron scheme. The k- 
point mesh and the plane wave kinetic energy cutoff (700 
eV) are chosen such, that total energies are converged 
to a numerical accuracy of 1 meV per formula unit.^'' 
The structures of UAUii,^^ NaAlH4,2i Mg(AlH4)2,^^ 
Ca(AlH4)2r^ Ca(BH4)2r^ KA1H4^ and of the alkah 
boranatesi^i are taken from the literature. We additionally 
relaxed the atomic positions, but the relaxations were 
small and had only a minor effect on the total energies. 
Our calculated values compare well to those obtained in 




Model (eV / AH4) 

FIG. 1: Model formation energies (eV/AH4) compared to 
DFT and experimental values.— 



previous calculations Details on Mg(BH4)2 will 
be published elsewhere.— 

Fig.[l]shows that, despite its simplicity, the model gives 
FEs that are in quantitative agreement with both the 
experimental, and the calculated DFT values. The rms 
deviation of the model with the experimental and the 
DFT values is 0.27 and 0.33 eV/AH4, respectively. Note 
that these numbers are comparable to the rms deviation 
between the experimental and the first principles values, 
0.19 eV/AH4, which represents the state-of-the-art. An 
obvious source of error is our neglect of the details of the 
crystal structure, e.g., by using an average Madelung con- 
stant in Eq. Changing the Madelung constant by 5% 
changes the FE of the alkali compounds by 0.4 eV/AH4 
and that of the alkaline earth compounds by 0.7 eV/AH4. 
As these numbers are larger than the rms deviation of the 
model, one can conclude that the details of the crystal 
structure are relatively unimportant. 

The model also seems to work reasonably well for 
some other boranates. The model FE for Sc(BH4)3 is 
0.37 eV/BH4 higher than the DFT value calulated by 
Nakamori et. al,^ which is within the rms error bar given 
above. The model FE for Zn(BH4)2 and CUBH4 are 0.6 
and 0.8 eV/BH4 higher than the DFT values, respec- 
tively. For cations with a nominal charge Zc — 4, such 
as Zr or Hf, the model breaks down. The model FE then 
deviates by 2.5 eV/BH4 from the DFT valuesi^ As can be 
seen from Eq. ([2]), for large Zc i?crys becomes sensitive 
to small changes in the Madelung constant and the ionic 
radius of the cation, or in other words, to the details of 
the crystal structure. 

Fig. [1] shows that boranates are generally more stable 
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FIG. 2: The contributions to the formation energy (eV/AH4) 
according to Eq. ((T|. 



than alanatcs. The origin of this stabihty can be analyzed 
by decomposing the FE into the contributions according 
to Eq. ([T|), which is shown Fig. [21 The differences in for- 
mation energy of the elements E'oicm are to a large degree 
compensated by the differences in the formation energies 
-Eions of the ions from the atoms. The ionic crystal energy 
Ecrys of the boranates is however significantly larger than 
that of the alanates, which results in a larger stability of 
the latter. This is a size effect since the BHj anions are 
significantly smaller than the AlH^ anions. 

It has been observed that the dissociation energies of 



complex alkali hydrides into simple alkali hydrides in- 
crease with the atomic number of the alkali atom.^^ For 
the FEs from the elements the overall trend is not that 
clear. i?ciom and -Eions both decrease with increasing 
atomic number, see Fig. [2l which increases the stabil- 
ity. However, this is almost compensated by i?crys, which 
increases with the cation radius re- 
in the alkaline earth series the FE decreases with the 
atomic number. The dominant effect is a decreasing 
£'ions, which is due to a decrease in the ionization po- 
tentials of the cations, see Table HI 

To summarize we constructed a model for the for- 
mation energies (FEs) of alkali and alkaline earth 
alanates and boranates from the elemental solids and II2 
molecules. The model is based upon ionic bonding be- 
tween metal cations and AIH4 or BH^ anions. It can 
be constructed using simple energy values that are avail- 
able in the literature and it does not make use of explicit 
crystal structure information. Compared to experimental 
values, the model FEs have a similar accuracy as calcu- 
lated DFT values. The trends in the FEs over the series 
of compounds can be analyzed in terms of the individual 
contributions to the model. 

The authors wish to thank R. A. de Groot for useful 
and stimulating discussions. This work is part of the re- 
search programs of 'Advanced Chemical Technologies for 
Sustainability (ACTS)' and the 'Stichting voor Funda- 
menteel Onderzoek der Materie FOM)', both financially 
supported by the 'Nederlandse Organisatie voor Weten- 
schappelijk Onderzoek (NWO)'. 



^ F. Schiith, B. Bogdanovic, and M. Felderhoff, Chem. Com- 

mun. 20, 2249 (2004). 
^ A. Ziittel, P. Wenger, P. Sudan, P. Mauron, and S. I. 

Orimo, Mater. Sci. Eng. B-Solid State Mater. Adv. Tech- 

nol. 108, 9 (2004). 
^ Y. Nakamori et al., Phys. Rev. B 74, 045126 (2006). 

* S. C. Chung and H. Morioka, J. Alloys Compd. 372, 92 
(2004). 

^ P. Vajeeston, P. Ravindran, R. Vidya, H. Fjellvag, and 
A. Kjekshus, Cryst. Growth Des. 4, 471 (2004). 
O. M. L0vvik, O. Swang, and S. M. Opalka, J. Mater. Res. 
20, 3199 (2005). 

P. Vajeeston, P. Ravindran, A. Kjekshus, and H. Fjellvag, 
J. Alloy. Compd. 387, 97 (2005). 

* M. Born, Verhandl. Deut. Physik. Ges. 21, 679 (1919). 

^ R. C. Weast (ed.) and M. J. Astle (ed.), CRC Handbook of 
Chemistry and Physics, (CRC press, New York, 62 edition, 
1982), F 190-201. 

^° M.P. Tosi, F. Seitz (ed.) and D. Turnbull (ed.), Sohd 
State Physics 16, (Academic press. New York and Lon- 
don, 1965), Chap. 1, 1-113. 
R. Shannon, Acta. Cryst. A 32, 751 (1976). 
M. Smith and G. Bass Jr., J. Chem. Eng. Data 8, 342 
(1963). 

P. E. Blochl, Phys. Rev. B 50, 17953 (1994). 
^'^ J. P. Perdew and J. A. Chevary and S. H. Vosko and K. 



A. Jackson and M. R. Pederson and D. J. Singh and C. 

Fiolhais, Phys. Rev. B, 46, 6671 (1992). 

G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 
" G. Kresse and J. Hafner, Phys. Rev. B 49, 14251 (1994). 
" G. Kresse and J. Hafner, Phys. Rev. B 47, 558 (1993). 

G. Kresse and J. Furthmiiller, Comput. Mater. Sci. 6, 15 

(1996). 

G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 
(1996). 

In both the DFT calculations and the model no zero point 
energies are taken into account. 

M. J. van Setten, G. A. de Wijs, V. A. Popa, and 
G. Brocks, Phys. Rev. B 74, in press, cont-math/0609189| 
(2006). 

^2 O. M. L0vvik, Phys. Rev. B 71, 144111 (2005). 

2^ K. Miwa et al., Phys. Rev. B 74, 155122 (2006). 

P. Vajeeston, P. Ravindran, A. Kjekshus, and H. Fjellvag, 
J. Alloy. Compd. 363, L7 (2004). 

25 A. Peles, J. A. Alford, Zhu Ma, Li Yang, and M. Y. Chou, 
Phys. Rev. B 70, 165105 (2004). 

2® M. J. van Setten, G. A. de Wijs, and G. Brocks, in prepa- 
ration (2006). 

^"^ M. E. Arroyo y de Dompablo and G. Ceder, J. Alloys 
Compd. 364, 6 (2003). 



